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1. Introduction 

Pion electromagnetic form factor Fy(q 2 ) is one of the fundamental observables in hadron 
physics. An analysis of experimental data based on chiral perturbation theory (ChPT) at two loops 
leads to a precise estimate of the charge radius (r 2 )y A detailed comparison of (r 2 )v between 
ChPT and non-perturbative calculations on the lattice may provide a good testing ground for recent 
lattice simulations in the chiral regime as well as a better understanding of the chiral behavior of 
Fv(q 2 ). 

While there is no experimental processes directly related to the scalar form factor F$(q 2 ), the 
chiral behavior of the scalar radius (r 2 )s is interesting, as it provides a determination of the LEC I4 
and has a 6 times enhanced chiral logarithm compared to (r 2 )v- A non-perturbative determination 
on the lattice is challenging, because we need to evaluate disconnected three-point functions. 

In this article, we update our analysis of F v (q 2 ) reported at the last conference [^] with doubled 
statistics, and present newly obtained results for Fs(q 2 ). These quantities are measured on gauge 
configurations of two-flavor QCD on a 16 3 x32 lattice generated with the overlap quark action along 
the fixed topology strategy The lattice spacing determined from the Sommer scale ro = 0.49 fm 
is a = 0.1184(21) fm. We refer the reader to Refs.[Q, |5|] for detailed setup and overviews of our 
production simulations. 

2. Measurement of pion correlation functions 

We measure pion correlators through all-to-all quark propagators Contributions of 100 
low-lying modes (A^, u^) (k= 1, ...,N ep ;N ep = 100) of the overlap operator D are evaluated ex- 
actly, whereas the remaining high modes are taken into account stochastically by the Z2 noise 
method. We prepare a single noise vector for each configuration, and dilute [||] it into Nj = 
3 x 4 x N t /2 vectors 77 W (k = l,...,Nd) with support on a single value for color and spinor in- 
dices and at two time-slices. The all-to-all propagator can then be expressed in a simple form 
D 1 = v(k) wWt (Nvec = N ep + N d ) with two set of vectors 

v« = i^,...,^,xW,...,x^, W ^ = {u^\...,u^\^\...,^}, (2.1) 

where =ZT 1 (1 u<M) r] W. 

From the v and w vectors, we may construct meson fields at a temporal coordinate t with the 
Dirac matrix T and spatial momentum p 



t;p) = £0(r)w(x + r,O Wt rv(x,O (/) e-' px . (2.2) 



For the smearing function 0(r), we choose the local 0/(r) = 8 r o and exponential function s (r) = 
exp[— 0.4jr|]. Connected and disconnected three-point functions as well as the subtraction term of 
the vev contribution to the scalar form factor, shown in Fig. [IJ, are calculated from these meson 
fields as 



C£?(*,*W) = i£ £ ^ + A, + A/;pOO + Ar;p-p') 

* t fe,/,m=l 

x^^-p), (2.3) 
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Figure 1: Connected (left-most diagram) and disconnected (middle diagram) three point functions. Note 
that we have the contribution to Fs(0) from the right-most diagram due to the non-zero vacuum expectation 
value of the scalar operator S. 



Nvec 

x^7»(HA(;p-p'), (2.4) 

m=l 

c£J(Af,Af';p,p') = l£ £ ^(t + At + At'- P ')^{t-- P ) 

*(^LL<rV;P-P')\ , (2-5) 

where (• • • ) con f represents a Monte Carlo average. We denote the temporal separation and spatial 
momentum for the initial (final) meson by At and p (At' and p'), respectively. 

Our measurements are carried out at four values of the quark mass m u d in the range 290 < 
Mx [MeV] < 520. We explore the region of the momentum transfer — l-l<q 2 [GeV 2 ] < by taking 
the meson momentum p with |p| <2. Note that the spatial meson momentum is shown in units of 
2na/L in this article. While we have simulated only the trivial topological sector, the effect of the 
fixed global topology is suppressed by the inverse of the space-time volume ~ 1/V [ff]. 



3. Determination of pion form factors 

We calculate effective value of the vector form factor from a ratio 

F (At At' ■ n 2 \ = 1Mn Ry(At,At';\p\,\p'W) 

W ,?J ^(|p|)+£*(|p'|) R v (At,At';0,0,0) ' 1 } 

* fAtAS\ \\ '\ ^ C^\At,At';p,p') 
R v (At,At;\p\,\p \,q z ) = , TK^ZTy 

Here C %% $y is the pion two-point function with the smearing function <p (<p') for the source (sink) 



operator, and it can also be calculated from the meson field Eq. ( ]2.2[ ). We take the average of Ry 
over momentum configurations corresponding to the same value of q 2 . This average as well as 



that over the location of the source operator in Eqs. ( p.3| ) - ( p.5[ ) leads to an accurate estimate of 
Fv(At,At';q 2 ) as shown in Fig. |2[ The vector form factor F n (q 2 ) is determined from a constant fit 
in a range of (At, At'), where Fy(At, At';q 2 ) shows a reasonable plateau. We include the leading 
finite volume correction (FVC) |Q] to Fy(q 2 ). 

The scalar form factor normalized at a certain momentum transfer q 2 ef can be calculated from 
a similar ratio 
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Figure 2: Effective value of F v (At,At';q 2 ) (left panels) and F s (At,At';q 2 )/Fs(At,At';q 2 ef ) (right panels) at 
' M wf/^' M s,phys/2, where m SiP h ys is the physical strange quark mass. 



F s (At,At';q 2 ) _ R s (At,At';q 2 ) 



R s (At,At';q 2 



C n i n (At,At';-p,p') 



(3.3) 



F s (At,At';q 2 ei ) /? s (Af,Af';«7 2 ef ) ' ' " , ' C Sff)0I0 ,(A/;p)C S3t , 0(0J (A/';p') ' 

where C^i^ =C^ 1 °™' ) -C^ ' +C^\ As Fig. | indicates, F 5 (^ 2 ) at q 2 =0 suffers from a relatively 
large statistical error than those at q 2 ^ due to the severe cancellation between cj^t and cj^j' . 
We therefore use Fs{q 2 ) normalized at the smallest non-zero momentum transfer with |q re f | = 1 in 
the following analysis. The normalized form factor Fs(q 2 )/Fs(q 2 ef ) is determined by a constant fit, 
while FVC to Fs(q 2 ) is not available so far and is not taken into account. 



4. Parametrization of q 2 dependence 

The vector and scalar form factors are plotted as a function of q 2 in Fig. ||. We observe that 
Fv{q 2 ) is close to the pole dependence 1/(1 —q 2 /M 2 ) with M p measured at simulated m uc j. Its 
q 2 dependence is therefore parametrized by the following form of the p pole with a polynomial 
correction to determine the charge radius (r 2 )v and the curvature cy 

Fv(q 2 ) = l _ q \ /M 2 + c ^ 2+C2 ^ 2 ) 2 + ^ (? 2 ) 3 = l + \{r 2 )vq 2 + c V {q 2 ) 2 + --- . (4.1) 

Because the deviation from the p pole is small, we obtain a reasonable % 2 /dot ~ 1, and results for 
(r 2 )v and cy are stable against the inclusion of the cubic correction term. 




« i 
Figure 3: Vector (left panel) and normalized scalar form factors (right panel) at m uc i^m s p h ys /2 as a function 
of q 2 . Solid lines show our parametrization and its statistical error. In the left panel, we also plot p pole 
contribution expected from the vector meson dominance hypothesis by the dashed line. 
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Figure 4: Chiral extrapolation of (r 2 )v (left panel) and (r 2 )s (right panel) based on one-loop ChPT. Star 
symbols show the experimental value for (r 2 )v [0] and an indirect determination of (r 2 }s through %% scat- 
tering rfi0|. 



Such pole contribution in the scalar channel is not clear within our statistical accuracy. Our 
data can be fitted to a simple quadratic form 

F s {q 2 ) = F S (0) (l + \(r 2 ) s q 2 +c s {q 2 ) 2 ^j . (4.2) 

We confirm that the result for the scalar radius (r 2 )s is stable if we switch to the cubic or single 
pole form F$(0)/(1 — g 2 /Mj? t ) with Mfi t as a fit parameter. The curvature cs is however strongly 
depends on the choice of the parametrization form, and hence is not used in the following analysis. 



5. Chiral extrapolation 



In one-loop ChPT, the radii (r 2 )y and (r 2 )s are given by 



(A 



1 



1 



NF 



NF 2 " b/ NF 2 

1/13 



In 



NF 2 



In 



M: 



21 



A 4 



(5.1) 
(5.2) 



where N = (47r) 2 . We set the renormalization scale jit to 4nF, and fix F to the value determined 
from our study of the pion mass and decay constant |§J. These fits are however not quite suc- 
cessful as seen in Fig. ||. While our data of {r 2 )v are fitted with reasonable # 2 /dof ~0.3, the 
value extrapolated to the physical quark mass 0.362(4) fm 2 is significantly smaller than experi- 
ment 0.437(16) fm 2 On the other hand, the one-loop formula for (r 2 )s with the enhanced 
chiral log fails to reproduce our data and results in large # 2 /dof ~ 16. We note that similar mild 
quark mass dependence of the radii is also observed by the ETM Collaboration with a different 



discretization on a slightly finer lattice []1 1[]. It is unlikely that the failure of the fits within one-loop 
ChPT is caused by systematic uncertainties due to the fixed topology and the finite lattice spacing. 
We then extend our analysis to two loops. The higher order contributions to the radii are given 

by§Bl 



A(r 2 ) y 



1 



N 2 F 4 \ 192 



13N 181 



+ 6^! 



19 



12M[ 2 ) M 2 ln 



Ml 



(5.3) 
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Figure 5: Simultaneous chiral extrapolation of (r 2 )v (left panel) and cy (right panel) based on two-loop 
ChPT (solid lines). We also plot one- and two-loop contributions by dashed and dotted lines, respectively. 
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At two loops, the curvature cy has non-trivial contributions and can be included in our analysis 
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(5.5) 



The analytic terms containing r r Xi (X = V,S, i= 1,2) represent contributions from 0(p 6 ) chiral 
Lagrangian. We denote the linear combination l\ — l r 2 j2 appearing commonly in (r 2 )v and cy as 



F 

'l,2- 

While the two-loop formulae involve many 
LECs, the simultaneous fit to {r 2 )v and cy has 
only four free parameters 1' 6 , l\ 2 , r r v j and r r V2 . This 
fit plotted in Fig. || shows that the two-loop contri- 
butions are significant in our simulated region of 
mud- We obtain a reasonable value of # 2 /dof~0.5, 
and the extrapolated values of (r 2 )y and cy are 
consistent with experiment [jl]|. 

The inclusion of {r 2 )s into the simultaneous 
chiral fit introduces additional four free parame- 
ters. In order to stabilize this fit, we fix V 2 and l^, 
which appear only in the two-loop terms, to a phe- 
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Figure 6: Chiral exttapolation of (r 2 )s from si- 
multaneous to (r 2 )v,s and cy based on two-loop 
formulae. 



nomenological estimate Z2 =4.31 [|10J and a lattice 
result Z3 = 3.44 from our analysis of the pion spectroscopy 



l . The extrapolation of (r 2 )y and cy 



turns out to be consistent with those in Fig. |5[ The extrapolation of (r 2 )s is shown in Fig.^. From 
this simultaneous fit, we obtain 



(r 2 )v 



: 0.404 (22) (22) fm 2 , (r 2 ) s = 0.578(69)(46) fm 2 , cy = 3. 1 1(14) (86) GeV -4 . (5.6) 



'The )i independent convention // is defined from /[ = )f(/; + ln[M|//x 2 ])/2A' with 73 = -1/2, 74 = 2 and 76 = — 1/3. 
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The first error is statistical. The second is systematic error estimated by changing the inputs for 
1' 2 and l\ to different phenomenological estimates in Ref. [[[]], and by limiting the fitting data to the 
radii ((r 2 )y and (r 2 )s) or those in the vector channel ((r 2 )y and cy). We also test (r 2 }s from the 
cubic parametrization for the q 2 dependence of Fs(q 2 ). Note that all the extrapolated values in 
Eq. ( ^6| ) are consistent with experiment. 

We obtain l 6 = 1 1 .8(0.7)( 1 .3), l 4 = 4.06(44) (99), and l\ 2 = -2.9(0.8)(2.4) x 10~ 3 for the 
0(p 4 ) LECs. Our estimate of l(, is slightly smaller than 16.0(0.9) obtained in Ref.[[l]] partly due 
to a deviation of F between our lattice determination ^ and two-loop ChPT [13]. We note that 
I4 is consistent with our determination I4 = 4. 12(56) from F % ^ and a phenomenological estimate 
4.39(22) 0]. Our results for the 0(p 6 ) LECs are r r Vi = -1.1 x 10~ 5 , r r v2 = -4.0 x 10~ 5 and 
rj j = 1.3 x 10~ 4 with substantial uncertainty of 50-100 %. 



6. Conclusions 



In this article, we report on our calculation of the pion form factors with two flavors of dynam- 
ical overlap quarks. By employing the all-to-all quark propagators, F${q 2 ) is calculated including 
contributions from the disconnected diagrams for the first time. The one-loop ChPT formulae fail 
to reproduce our data of Fs(q 2 ). In our analysis extended to two loops, we observe significant 
two-loop contributions at our simulated quark masses, and obtain (r 2 )ys and cy consistent with 
experiment. Further investigations of systematics due to the fixed global topology and quenching 
of strange quarks are in progress by direct simulations in the non-trivial topological sectors and in 
three-flavor QCD. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 08-05). This work is supported in part by the Grant-in- Aid of the Ministry 
of Education (No. 18340075, 18740167, 19540286, 19740160, 20025010, 20039005, 20340047, 
and 20740156), the National Science Council of Taiwan (No. NSC96-2112-M-002-020-MY3, 
NSC96-2112-M-001-017-MY3, NSC97-2119-M-002-001), and NTU-CQSE (No. 97R0066-69). 
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